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Free expansion of attractive and repulsive Bose-Einstein condensed vortex states formed in an axially 
symmetric trap is investigated using the numerical solution of the time-dependent Gross-Pitaevskii 
equation. In a repulsive condensate the vortex-core radius is much smaller than the radial root- 
mean-square (rms) radius, which makes the experimental observation of the vortex core difficult. 
The opposite is found to be true in an attractive condensate which makes it a better candidate 
for experimental observation. Also, in all cases the ratio of vortex-core radius to radial rms radius 
increases as the angular momentum of the vortex increases. This makes the vortex states with 
higher angular momenta more suitable for experimental confirmation. 
PACS Number(s): 03.75.Fi 
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I. INTRODUCTION 

Since the successful detection [Q of Bose-Einstein con- 
densates (BEC) in dilute bosonic atoms employing mag- 
netic trap at ultra-low temperature, one problem of ex- 
treme interest is the formation of vortex states in an ax- 
ially symmetric trap. The recent experimental detection 
H of vortex states in the condensate has intensified this 
endeavor. There have been many theoretical studies on 
different aspects of BEC ||-[| and specially, on the vortex 
states in axially symmetric traps |q-[l3||. 

The study of the free expansion of a vortex state is 
of utmost importance as such an expansion can aid in 
its experimental confirmation [Q. The specific behavior 
of the time evolution of vortex core and the radius of 
the atomic core after the removal of the trap could be 
fundamental in the identification of vortex states. How- 
ever, a previous theoretical study on the free expansion 
of vortex states was limited to condensates with repul- 
sive interatomic interaction |Q. In this work we extend 
this effort to condensates with attractive interatomic in- 
teraction and compare critically the expansion dynamics 
with that for repulsive interaction. This is of extreme 
experimental relevance as by changing the external elec- 
tromagnetic field now it is possible to change the effective 
interatomic interaction in a controlled way from repulsive 
to attractive E3 |. 

We base our study on the numerical solution P,|o],[7|p3| 
of the nonlinear time-dependent mean-field Gross- 
Pitaevskii (GP) equation |H| with an axially symmetric 
trap, which should provide a faithful description of the 
formation and evolution of vortex states |7] . 

The study of superfluid properties of BEC is of great 
interest to both theoreticians and experimentalists . 
Quantized vortex state in BEC is intimately connected to 
the existence of superfluidity Such vortices are expected 
in superfluid Hejp However, due to strong interaction 
between the helium atoms there is no reliable mean-field 
description and a controlled theoretical study of the dy- 
namics is not possible. 

On the other hand, a weakly interacting trapped BEC 



is well-described by the mean-field GP equation that is 
known to admit vortex solutions for a trap with cylindri- 
cal symmetry which can be studied numerically. In 
contrast to superfluid Hejj, the condensed vortex state of 
trapped BEC is an extremely dilute quantum fluid and 
permits a mean-field description. Moreover, vortex states 
in BEC have been observed experimentally [Q . Many dif- 
ferent techniques for creating vortex states in BEC have 
been suggested ||, e.g., stirring the BEC with an exter- 
nal laser B , spontaneous vortex formation in evaporative 
cooling p| , and rotation of axially symmetric trap |l2[ . 

It is possible to have dynamically stable vortex BEC 
states with low quanta of rotational excitation or angu- 
lar momentum L per particle for both attractive 
and repulsive JH|] atomic interactions. However, in the 
attractive case such states exist for the number of atoms 
below a critical number which increases with L. In the 
absence of vortex (L — 0), the stable condensate in a 
cylindrical trap has a cylindrical shape. Such a BEC has 
the largest density on the axis of the trap. In the pres- 
ence of vortex motion the region of largest density of the 
BEC with nonzero L is pushed away from the central 
axial region and the atoms have more space to stabilize. 
The vortex state of the condensate in a cylindrical trap 
has the shape of a hollow cylinder with zero density on 
the axis of symmetry. 

For vortex states with low angular quanta L and large 
number of atoms, the radius of the hollow part is much 
smaller than the radius of the condensate in the radial 
direction, which makes the experimental confirmation of 
vortex states difficult. It has been argued and estab- 
lished in repulsive condensates with L = 1 that after a 
free expansion, the hollow core of the vortex may expand 
faster than the radius of the condensate and thus make 
the internal hollow region detectable experimentally con- 
firming a vortex state . 

In our study of free expansion of vortex states we find 
that the attractive condensates allow vortex states with a 
different nature from those in repulsive condensates. For 
a fixed L, the ratio of vortex core radius to radial root- 
mean-square (rms) radius of an attractive condensate is 
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always larger than the same ratio for a repulsive con- 
densate, which may make attractive condensates more 
appealing for generating, detecting and studying vortex 
states. We also find that for a fixed L, upon free expan- 
sion the above ratio increases for repulsive condensates, 
whereas it decreases for attractive condensates remain- 
ing, however, always larger than the corresponding ratio 
for repulsive condensates. In all cases the ratio increases 
with L. 

The present study is performed with the direct numer- 
ical solution of the time-dependent GP equation with an 
axially symmetric trap. In the time-evolution of the GP 
equation the radial and axial variables are dealt with in 
two independent steps. In each step the GP equation 
is solved by discretization with the Crank-Nicholson rule 
complimented by the known boundary conditions Jl3| , ^6| . 
We find that this time-dependent approach leads to good 
convergence. 

In Sec. II we describe briefly the time-dependent GP 
equation with vortex states and a numerical method for 
its solution. In In Sec. Ill we report the numerical results 
of the present investigation about the free expansion of 
vortex states and finally, in Sec. IV we give a summary 
of our investigation. 



II. NONLINEAR GROSS-PITAEVSKII 
EQUATION 

The deatil of the calculational procedure based on the 
GP equation has been elaborated elsewhere jl3| and we 
present a summary of our account here. At zero tem- 
perature, the time-dependent Bose-Einstein condensate 
wave function ^(r;r) at position r and time r may be 
described by the following mean-field nonlinear GP cqua- 
tion |,[l§ 
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*(r;r) = 0. 



(2.1) 



Here m is the mass and N the number of atoms in the 
condensate, g — Airh 2 'a/m the strength of interatomic in- 
teraction, with a the atomic scattering length. The trap 
potential with cylindrical symmetry may be written as 
V(r) = ^muj 2 (r 2 + X 2 z 2 ) where u> is the angular frequency 
in the radial direction r and \u> that in the axial direc- 
tion z. We are using the cylindrical coordinate system 
r = (r, 0, z) with 9 the azimuthal angle. The normaliza- 
tion condition of the wave function is J dr|^(r; t)\ 2 = 1. 

The GP equation (2.1) can accommodate quantized 
vortex states with rotational motion of the condensate 
around the z axis. In such a vortex the atoms flow 
with tangential velocity Lh/(mr) such that each atom 
has quantized angular momentum LTi along z axis. This 
corresponds to the angular dependence 



of the 



^(r; r) = tp(r, z; r) exp(iL#) 
\ function. 



(2.2) 



Substituting Eq. (2.2) into Eq. (2.1), and transform- 

V2r/l, 



ing to dimensionless variables defined by x 
y = \/2z/l, t = tuj, I = yjfij (muj), and 
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<fi(x,y;t)=0, (2.4) 



where n = Na/l. The normalization condition of the 
wave function become 
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For a stationary solution the time dependence of the 
wave function is given by (p(x,y;t) — cxp(—i[it)ip(x,y) 
where /i is the chemical potential of the condensate in 
units of huj. If we use this form of the wave function in 



Eq. (2.4), we obtain the following stationary nonlinear 



time-independent GP equation |15j 
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<p(x,y) = 0. (2.6) 



The rms radii of the condensate in the radial and axial 
directions are defined, respectively, by 

/>oo />oo 

(zrms) 2 = 2tt / dx dyx\tp{x,y;t)\ 2 , (2.7) 
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and 



(yrms) 2 = 27r/ dx dyx^y 2 \<p(x, y; t)\ 2 . (2. 



We solve the GP equation (|2.4|) using a time-iteration 
method elaborated in Refs. [|1 3j|l 6f| - The full GP Hamil- 
tonian is conveniently broken into two parts H x and H y 
— the first containing the s-dependent terms and the sec- 
ond containing the y-dependent terms with the nonlin- 
ear interaction term divided equally into two parts. The 
GP equations along the x and y directions are defined 
on a two-dimensional set of grid points N x x N y using 
the Crank-Nicholson discretization method. The resul- 
tant tridiagonal equations along x and y directions are 
solved alternately by the Gaussian elimination method 
along the x and y directions [M. Effectively, each time 
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iteration of the GP equation is broken up into two parts 
— first using H x and then H y . For a small time step A 
the error involved in this break-up procedure along x and 
y directions is quadratic in A and hence can be neglected. 
This scheme is repeated for about 500 time iterations to 
yield the final solution. 

The initial solution for the iterative scheme is taken 
to be the analytic solution of the corresponding free har- 
monic oscillat or s etting the nonlinear term to zero in the 
GP equation ( |2.6| ). At each iteration the strength of the 
nonlinear term is increased by a small amount so that af- 
ter the final iteration the full strength is attained and the 
required solution of the GP equation obtained. The free 
expansion of the vortex states is studied essentially using 
the above iterative scheme by setting the trapping po- 
tential to zero in the GP equation after the final solution 
for the vortex state is obtained. 

The typical value of the space step used for discretiza- 
tion along x and y directions is 0.02 and that for time 
step A is 0.05. For small nonlincarity the largest values 
of x and y are x max = 8, |y|max = 8. However, for 
stronger nonlinearity, larger values of £max and |y|max 
(up to 20) are employed. 

III. NUMERICAL RESULT 

We present results for the free expansion of vortex 
states using the numerical solution of the time-dependent 
GP equation for attractive and repulsive interatomic in- 
teractions for different L. In a recent study of repulsive 
vortex states with L = 1, Dalfovo et al. pointed out 
that the small vortex core radius of these condensates 
make the vortex states hard to detect experimentally. 
They show that after a free expansion, the vortex-core 
radius increases faster than the radial rms radius, which 
makes the experimental observation of the vortex core 
easier. Here we complement the above study by includ- 
ing the attractive condensates as well as vortex states 
with L > 1. The interesting features of the results that 
we emphasize here are quite independent of the value of 
axial to radial trap frequency A, and in this study we 
present results for A = \/8 only. 

In Figs. 1(a) and (b) we plot the profile of the section 
of the wave function | <^>(ar, 0; t) | = \ip(x,0;t)\/x upon free 
expansion as a function of x at different times < t < 3 
for n — —1.5 (attractive case) and n = 12 (repulsive 
case), respectively, — the trap being removed at time 
t = 0. In the attractive case, the minimum value of 
n allowed for L = 1 is —1.63 A vortex-core ra- 

dius x c is conveniently defined to be the smaller value 
of x where the central core density \<fi(x, 0; t)\ 2 attains 
er l times the peak value in Figs. 1(a) and (b). After 
free expansion, both the core radius x c and radial rms 
radius of the condensate increase. The relative increase 
of these two radii is crucial in the experimental observa- 
tion of the vortex states. For this a careful study of the 



evolution of the different radii is important. In Fig. 2 
we present the time evolution of irms and yrmsj where 
we plot xrms(i)/^rms(0) and j/rms(*)/yrms(0) vs. t for 
L = 1 and n = —1.5,0,12. Both ratios increase with 
time in all cases. In the attractive case the radius in the 
radial direction increases faster, the opposite being true 
in the repulsive case. The result for n = always stays 
between the attractive and repulsive cases. 

Next we study the evolution of the ratio ic/irmsj al- 
ways less than unity, with time t. The larger this ratio, 
the larger is the vortex-core radius compared to the radial 
rms radius and the easier could be the experimental de- 
tection of the vortex state. In Fig. 3 we plot x c /xrms vs. 
time t for L = 1 and 4 for several n including the attrac- 
tive (negative n), repulsive (positive n) and noninteract- 
ing (n = 0) cases. For fixed n and t, this ratio increases 
with L; and for fixed L and t, the ratio increases as n 
decreases from positive to negative values (passes from 
repulsive to attractive). This means that for attractive 
interaction the vortex state has a larger vortex-core ra- 
dius compared to the radial rms radius and hence could 
be easily detected. From Fig. 3 we find that at t = 
the numerical value of this ratio for attractive condensate 
could be twice as large as its value for repulsive conden- 
sate. However, there is an interesting difference between 
the time evolution of this ratio for a fixed L for attrac- 
tive and repulsive interactions. This ratio decreases with 
time for attractive interaction and increases with time for 
repulsive interaction. Hence there is no specific advan- 
tage in allowing an attractive vortex to expand for the 
sake of detection. It will be easier to detect it before ex- 
pansion, when the vortex core is more pronounced. For 
the noninteracting case the above ratio remains constant 

at all times 0.28 for L = 1 and - 0.60 for L = 4. 

Hence, although the possibility of detection of the repul- 
sive vortex state increases upon free expansion, the op- 
posite is true for attractive interaction. These features of 
the attractive vortex states upon expansion are the most 
important results of this investigation clearly displayed 
in Fig. 3. 

We also studied the expansion of vortex states of higher 
angular momenta in detail. The general features of these 
states are similar and we present the results for L = 4 
in the following. In Figs. 4 (a) — (d) we plot the profile 
of the complete wave function \4>{x, y; t)\ vs. x and y for 
n = —4.3 for t = 0,1,2 and 3. In this attractive case the 
critical n for collapse is —4.4 In Figs. 5 and 6 we plot 
the same for n = 4.3 and 16, respectively. At all times the 
wave function is zero on the y axis and is peaked at some 
finite x. Figures 4 — 6 clearly show the time evolution 
of the full condensate in these cases. We find from these 
figures that at a particular time, the vortex core radius x c 
is approximately the same in the three cases. However, 
the radial rms radius x-^ms increases as n increases from 
negative to positive values and consequently, the ratio 
ic/irms increases as we pass from Fig. 4 to 5 and then 
to 6 (from attractive to repulsive). The evolution of the 
ratio x c /xrms upon free expansion for L — 4 and different 
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n is shown in Fig. 3. The results for the L = 4 case is 
qualitatively similar to the L — 1 case and the discussion 
and comments of the L = 1 case apply here. 

An interesting feature of the freely expanding conden- 
sate emerges from the present study. The vortex core 
radius x c and the radial rms radius Xyms play important 
roles in this game. The condensed vortex state has the 
form of a hollow cylinder. Qualitatively, xrms is the ra- 
dius of the cylinder and x c is the radius of the hollow. 
The hollow space is more prominent when x c /xrms is 
larger. A large value of x c /x r ms facilitates the experi- 
mental confirmation of the vortex state. From Figs. 1, 
2, 4, 5 and 6 we find that both x c and irms increase 
upon free expansion. However, for an attractive conden- 
sate srms increases at a faster rate than x c so that the 
above ratio decreases with time. The opposite is true for 
a repulsive condensate. For a noninteracting condensate 
both radii increase at the same rate and the ratio remains 
constant. The whole scenario is physically understand- 
able. For an attractive condensed vortex state upon free 
expansion, because of interatomic attraction, the radius 
of the hollow x c can not increase fast enough with time 
compared to the radius of the cylinder. The cylinder as a 
whole swells at a faster rate than the hollow making the 
ratio cCe/xrms decrease with time. The opposite happens 
for a repulsive condensate, where the atomic repulsion 
makes the hollow part swell at a much faster rate. For 
a noninteracting condensate x c and Xrms increase at the 
same rate leading to a constant x c /xrms- In general for 
a fixed L, a; mis is much smaller for an attractive conden- 
sate, which makes the ratio x c / x^ms much larger than 
the same for a repulsive condensate. 



IV. SUMMARY 

In this paper we present a numerical study of the time- 
dependent Gross-Pitaevskii equation for an axially sym- 
metric trap to obtain insight into the free expansion of 
vortex states of BEC. The time-dependent GP equation 
is solved iteratively by discretization using a two-step 
Crank-Nicholson scheme 
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The ratio of vortex core radius to radial rms radius 
%c/ Xrms(< 1) is found to play an interesting role in the 
free expansion of condensed vortex states. The larger 
this ratio, more prominent is the vortex core and the 
more easier is the possibility of experimental detection of 
vortex states. For stable vortex states of fixed angular 
momentum per particle L, this ratio for attractive inter- 
action is larger than that for repulsive interaction. Upon 
free expansion of such states, the ratio for attractive 
interaction decreases remaining, however, always larger 
than that for repulsive interaction, which increases after 
expanion. The ratio increases with L in all cases. This 
makes the stable attractive vortex states before expan- 
sion more suitable for experimental detection. 
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Figure Caption: 

1. Wave function 4>(x,y — 0) of a stable vortex state 
with L = 1 vs. x for (a) n = —1.5, and (b) n = 12. 

2. The ratios x rms (t) / x rms (0) and 2/rms(i)/j/rms(0) 
vs. time t of expanding vortex states with L = 1 and 
n = -1.5,0,12. 

3. The ratio x c (t)/xrms{t) vs. time t of expanding 
vortex states with L = 1, n = —1.5,0,12 and L = 4, 
n = ±4.3,0,16. 

4. The wave functions 4>(x, y) of a freely expanding 
vortex state vs. x and y at times t = 0, 1, 2, 3 for L = 4 
and n = —4.3. 

5. Same as Fig. 4 for n = 4.3. 

6. Same as Fig. 4 for n = 16. 
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